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ABSTRACT 


A  Fortran  computer  program  has  been  developed  which  ana¬ 
lyzes  a  nonlinear  material  to  determine  (a)  the  range  of  IR  wave¬ 
lengths  that  can  be  converted  In  a  phase  matched  (PM)  process, 
(b)  the  PM  orientation  of  the  wave  vectors  for  critical  and  non- 
crltlcal  PM,  (c)  the  angular  aperture  for  PM  conversion,  and  (d) 
the  maximum  number  of  resolvable  lines  for  Image  conversion. 
These  characteristics  are  determined  as  a  function  of  IR  wave¬ 
length  for  a  given  pump  wavelei^^,  pump  radiation  divergence, 
and  length  of  nonlinear  crystal. 
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A  COMPUTER  PROGRAM  TO  ANAI YZE  OPTICAL  PARAMETRIC 
UP-CONVERSION  PROCESSES  IN  NONLINEAR  CRYSTALS 


INTRODUCTION 

Infrared  (IR)  light  can  be  up-converted  to  visible  light,  via  a  parametric  process,  in 
crystals  which  have  a  nonlinear  auscepUbillty  (1-3).  In  the  second-order  processes  con- 
si^red  here,  IR  light  of  frequency  Is  mixed  with  intense  light  of  frequency 
called  the  pump  light.  The  result  is  signal  light  of  frequency  For  the 

process  to  be  efficient,  the  crystal  must  be  transpai'ent  at  the  three  frequencies  involved, 
and  momentum  must  be  conserved;  l.e.,  the  process  must  be  "phase  matched"  (PM).  A 
process  Is  PM  if  the  directions  of  propagation  in  the  nonlinear  crystal  for  the  three  fre¬ 
quencies  of  light  can  be  found  such  that  the  signal  wave  vector  equals  the  sum  of  the  pump 
and  IR  wave  vectors,  l.e.,  (k^  +  ^jg)  -  =  0. 

Given  a  nonlinear  crystal.  It  Is  nscessary  to  find  the  range  of  IR  wavelengths  that 
can  be  up-converted  In  the  material  and  to  find  the  characteristics  of  the  up-converslon 
process  as  a  function  of  the  IR  wavelength.  Up-converslon  of  IR  to  the  visible  necessi¬ 
tates  the  use  of  an  intense  source  of  laser  light  as  a  pump  whose  wavelength  Is  near  the 
long-wave lengtti  edge  of  the  visible  spectrum.  Once  the  wavelength  of  the  pump  source 
has  been  specified,  the  range  of  ujg  for  PM  up-converslon  can  be  determined.  Then  for 
each  set  of  frequencies  tog,  ajg,  and  us,  the  orientation  of  the  crystal  with  respect  to 
the  wave  vectors  that  gives  a  PM  up-converslon,  l.e.,  the  phase  match  angles,  can  be 
found.  Jn  each  case  the  wave  vectors  can  be  either  colUnear  or  noncolUnear,  and  both 
sltmtlons  must  be  considered.  Also,  If  the  IR  to  be  up-converted  carries  spatial  Infor¬ 
mation  such  as  an  Image,  one  must  determine  such  things  as  the  acceptable  angular  ap¬ 
erture  for  PM  up-converslon,  the  resolution  limit*  for  the  beam  divergence  of  a  given 
pump,  and  the  orientation  which  maximizes  the  angular  aperture.  Further,  all  of  the 
characteristics  of  a  PM  up-converslon  process  must  be  calculated  as  a  function  of  ujg, 
the  IR  frequency. 

Finding  the  PM  angles  and  other  characteristics  of  an  up-converslon  process  means 
solving  transcendental  equations.  To  obtain  any  desired  accuracy  the  complete  analysis 
is  best  done  on  a  computer  using  an  Iteration  technique.  This  report  describes  a  Fortran 
rv  program  w!.lch  does  (hat  analysis.  The  characteristics  of  the  nonlinear  crystal  nec¬ 
essary  for  the  analysis  are  Us  limits  of  transparency  and  Its  refractive  Indexes  as  a 
function  c'  wavelen^. 

The  next  section  describes  the  mathematical  basis  of  the  calculation.  This  Is  fol¬ 
lowed  by  a  description  of  the  program  itself  and  two  appendixes  which  give  a  listing  of 
the  program,  sample  input  data,  and  the  corresponding  output. 


MATHEMATICAL  FORMULATION 

The  requirements  for  PM  up-converslon  (4)  are 

Vf  +  ujg  =  «(5 

and 


(1) 


*The  resolution  limit  considered  here  is  the  particular  case  of  an  extreme  multimode  pump  source 
with  the  IR  Image  and  pump  source  optiosUy  at  Infinity  with  respect  to  the  nonlinear  crystal.  See, 
for  example,  *IR  Image  Optical  Parametric  Up-ConverBion,”  R.A.  Andrews,  IEEE  J.  Qusnt.  Eleotr., 
Jen.  1970. 
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kr  +  k„  =  ks  .  (2) 

where  the  subscripts  F,  ig,  and  s  refer  to  pump,  Infrared,  and  signal,  respectively,  and 
k  Is  the  corresponding  wave  vector.  For  a  cryst^  of  length  L,  PM  up-converslon  takes 
place  as  long  as 

All  a  |Ak|  =  |k^  *  kjj,  -  kjl  S  .  (2) 

where  for  convenience  we  assume  that  Ak  is  measured  along  the  direction  of  ks  as 
shown  In  Fig,  1.  In  general,  the  orientation  of  the  nonlinear  crystal  for  PM  up-conversion 
is  given  by  the  solution  to  the  following  equation: 


Ak  = 


r>‘(up,$p)  in) 


sec  p  - 


(4) 


Where 


4  “  angle  between  kp  and  k^y; 

p  a  angle  between  kp  and  ks 

e,  a  angle  between  the  2-axls  and  kf<f  =  F,  ir,  S) 

i ,  > ,  k  a  0  or  £  for  extraordinary  or  ordinary  polarization 

It  Is  assumed  that  the  parametric  process  is  confined  to  the  yz  plane,  where  x,  y,  and 
t  are  the  optical  axes  of  the  nonlinear  crystal,  and  z  Is  the  optical  axis  for  a  uniaxial 
crystal.  However,  the  crystal  may  be  optically  biaxial;  therefore 


z 


Fig,  1  -  Crj’stal  geometry  for 
up-oonyerslon  from  infrared 
to  vlBible  wavelengths 


NRL  REPORT  6941 


3 


n°(<u,9)  =  n,(o) 


(5) 


and 


•in’  6  cot’  $ 
- + - - 

[ny(tt)))*  [n,(tu)]’ 


(8) 


Where  n^,  n^,  and  n,  are  the  refractive  Indexes  for  light  polarized  along  the  x,  y,  and 
I  directions. 

Solutions  to  Eq.  (4)  and,  in  particular,  colllnear  up*converslon  (^2  ^  s  o)  represent 
examples  of  critical  phase  matching.  The  noncrltlcal  phase  matching  (NCP^  solution 
(Ref.  r)  requires  that 


=  0  (7) 

The  crystal  has  finite  length  L  and  up-converslon  will  take  place  whenever  <  2k /L. 
Hence,  for  a  fixed  pump  beam  direction,  there  will  be  a  range  of  values  of  ^  for  which 
IR  light  will  be  up-converted.  This  range  defines  the  angular  aperture  s  for  and  Is 
defined  by  all  values  of  such  that  for  a  fixed  value  of  fff, 


<  2n/L  (8) 

or 

«  =  l«Si  -  <ii(  .  (9) 

where 


-  2ir/t  .  ^  ^0 

2  2f,/L  . 

2  0  . 


is  the  angular  aperture  Inside  the  crystal.  Outside  the  crystal, 

“fir  =  I  ^1)  »in  Sy  +  *a)  ^j)  ]  (10) 

In  this  manner,  each  PM  process  can  be  characterized  by  the  width  of  the  a*(0)  vs 
curve  with  dp  fixed.  Hence,  the  half- widths  are  defined  as 


and 


HW I  -  <t>^  ~  I 

HW ^  ~  . 


(11) 

(12) 


The  pump  radiation  always  has  a  finite  divergence,  and  hence  the  PM  angles  are  not 
well  defined.  Variations  In  Sp  will  give  PM  up-convorslon  for  cUlferent  values  of  0; 
hence  the  angular  apei'ture  Is  Increased.  The  dlverg<  nee  A  of  the  pump  radiation  Inside 
the  crystal  Is,  however,  less  than  that  outside,  l.e.. 
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^lyr  "  ■ 


(13) 


vher«  h  Is  some  average  Index  of  refraction. 


The  maximum  angular  aperture  Is  obtained  with  NCPM.  The  aperture  is  the  greatest 
U,  in  this  case,  Eqs.  (4)  and  (7)  are  solved  with  a*  =  -2n/L.  This  technique  Is  discussed 
by  Warner.  However,  a  larger  angular  aperture  does  not  always  guarantee  better  reso¬ 
lution.  To  determine  resolution,  one  must  find  the  change  In  <(,  that  will  cause  a  variation 
in  9s  Just  equal  to  the  width  of  values  caused  by  the  divergence  of  the  pump  radiation. 
A  varlfdion  in  Of  due  to  a  divergence  of  a  causes  a  change  In  i9s  of 


y  =  A  + 


sin  A) 


■ 

•  sin~* 

—  sir.  ,j> 

.  "s  J 

(14) 


This  change  In  65  is  equivalently  produced  by  a  change  In  (with  6p  fixed)  of 


( : —  sin 

A  ♦  sin**  (~  3 in 

V  h  J 

j 

(15) 


Hence,  the  number  of  resolvable  lines  of  IR  Is  equal  to 

^  ~  (1®) 


QENEIIAL  DESCRIPTION  OF  THE  PROGRAM 

The  program  is  written  In  Fortran  IV  language  and.  In  the  form  given  In  Appendix  A, 
hM  been  run  on  a  CDC  3800  computer.  Input  to  the  program  is  minimal  and  is  on 
punched  cards. 

..  or  each  material  to  be  analyzed,  it  Is  first  necessary  to  fit  dispersion  data  to  a 
Sellmeier  equation  of  the  form 

n/  =  .4^  +  B./(Cj -A^)  -  .  (i^x.y.z).  (17) 

(The  program  can  easily  be  modified  to  use  any  other  dispersion  relation  by  modifying 
the  function  FIN  (W,  I),  which  is  given  in  Table  1.)  Besides  the  constants  Ai,  b,,  c,  , 
and  Bj,  it  is  necessary  to  supply  the  upper  and  lower  wavelength  limits  of  the  transpar¬ 
ent  region  of  the  crystal,  the  length  of  the  crystal,  and  the  range  of  values  of  the  pump- 
source  divergence  to  be  considered.  This  Information  is  given  on  a  set  of  four  data 
cards  for  each  material.  The  required  format  Is  given  in  Table  2.  The  first  data  card 
gives  the  number  of  materials  to  be  analyzed.  Data  cards  2  through  5  are  repeated  for 
each  material.  The  symbols  used  lor  the  data  are  defined  in  Table  3. 

Pump  wavelengths  of  0.6943  m  and  1.06  u  have  been  selected,  since  they  He  in  or 
near  the  red  part  of  the  spectimm  as  discussed  above,  and  Intense  laser  sources  are 
readily  available  at  these  wavelengths.  The  program  analyzes  up-conversion  processes 
in  which  all  three  wavelengths  He  in  the  region  of  transparency.  The  IR  wavelength  Is 
scanned  from  the  lower  wavelength  Umit  of  transparency  upward  In  equal  increments  on 
a  Icgarithmic  scale.  The  Interval  Is  determined  in  statement  10  of  the  main  program. 
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Table  1 
Subprogram* 


Nam* 

Calculation* 

Lin*  No.  of 
Accuracy  Teat 

FNX  (XX,  PP,  DD) 

Angle  between  and  t-axl*  tor  PM 

where 

XX  -  angle  between  k^  and  k,, 

PP  ■  (*ee  text) 

DD  >  value  of  nk  at  PM  (eee  text) 

222 

TNCPM  (TCL,D) 

Noncritlcal  PM  angle 

TCL  -  coUlnear  PM  angle 

D  B  value  of  Ak  at  PM 

144 

PVAR  (T,  P,  DD,  B) 

Variation  in  ^  for  PM  when  9j>  1* 
varied 

T,P  -  PM  angle* 

DD  B  amount  is  increased 

E  B  value  of  Ak 

178 

FOK  (T,  P) 

Ak 

T,  P  s  PM  angle*  a^,  ^ 

- 

HWDK(L,T,  P,I) 

Width  of  Ak(^)  v*  ^  curve 

L  B  length  of  crystal 

T,  P  B  PM  angle* 

I  B  index  to  epecify  type  of 

variation  in  ^  to  be  performed 

281 

OUT 

Printout  routine 

— 

FIN  (W,I) 

Refractive  index 

W  B  wavelength  of  deeired  index 

I  -  1,  2,  or  3  correepondlng  to  x,  y, 
or  <  optical  axis 

■ 

1 

Refractive  index  for  arbitrary 
direction  of  propagation 

A  B  angle  between  direction  of 
propagation  and  t  axis 

[■ 

Refractive  index  sqiuared  for  arbi¬ 
trary  direction  of  propagation 

Q  -  eame  a*  A  in  FNl  (A) 
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Table  2 

Input  Data  Cards 


Data  Card 
Number 

Format 

Symbols 

1 

(Q) 

N 

2 

(AS,  5X,  6F  10.1) 

ANAME,  UL,  LL,  L,  DL,  DU,  DD 

3 

(4F  10.4) 

AX,  BX,  CX,  DX 

4 

(4F  10.4) 

AY,  BY,  CY,  DY 

5 

(4F  10.4) 

AZ,  BZ,  CZ,  DZ 

Table  3 

Definitions  of  Symbols 


Symbol 

Definition 

Units 

N 

Number  of  materials  to  be  aitalyzed 

Integer 

ANAME 

Name  of  material 

UL,  LL 

Upper  and  lower  wavelength  limits  for 
transparency 

Angstrom 

L 

Length  of  crystal 

Centimeter 

AX,  BX,  . . .  DZ 

Parameters  for  Sellmeler  equation  for  in¬ 
dexes  corresponding  to  x,  y,  and  z  opti¬ 
cal  axes.  (For  uniaxial  crystal  z  _  ex¬ 
traordinary  polarization,  x  =  y  —  ordinary 
polarization.) 

For  each  set  of  possible  wavelengths,  ail  polarization  combinations  are  checked  for  a 
possible  PM  process.  The  output  is  labeled  ”0  +  o  =  etc.  for  (ordinary  polariza¬ 
tion)  +  (ordinary  polarization)  =  (extraordinary  polarization),  etc.  "Not  Phase m ate hable" 
is  printed  if  an  orientation  cannot  be  found  for  phase-matched  up-conversion.  If  up- 
conversion  at  a  given  set  of  wavelengths  is  phase  matchable,  then  all  the  parameters 
discussed  in  the  previous  section  are  calculated  and  printed  (see  Appendix  B). 

The  functions  of  the  various  subprograms  are  listed  in  Table  1.  In  particular,  the 
function  FNX  (X,  P,  D)  is  more  general  than  need  be  for  this  program.  FNX  will  calcu¬ 
late  a  PM  angle  for  any  orientation  of  kji;  with  respect  to  k^.  In  this  case  P  is  the 
angle  between  the  (kp,  i-axis)  plane  and  the  (k,jf,  z-axis)  plane,  and  X  is  the  angle  be¬ 
tween  k/.  and  kjK- 

The  various  calculatio  ts  are  iterated  until  a  predetermined  accuracy  is  reached. 
The  line  numbers  for  the  appropriate  accuracy  determining  points  in  the  program  are 
listed  in  Table  1 . 
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appendix  A 


PROGRAM  LISTING 

PROGHAM  PhaSMCM  00000 lOO 

DATA  (UP= 5 7. 295/7951 .(  ON« ♦ 0 1 7453293 ),( CC-6. 283 1  868 »  » (Cl *1000 . 1  00000200 

COMMON/' 1 /W i  .W2tW3  /7/IN3/8/|N2/9/'|Nly'4/IX/10XPNC  00000300 

1  1  1  ,'K  .  TCU  •  TNCC  I  MWC  1  4  MWC2  •  HWNC  1  •  HWNC2  •  HWC3  •  HWC4  ■  MWNC3  .  HWNC4  .0  .  J  .0  I  V00000400 


2/1  2/CXHVii  1  •EXMK.'2 4EXm1<EX\ii2  •  EXHuS •EXHw6»FNL  1  ♦  TNLH  •PC  1  iPNC  1  •  PC2 

3  /5/MWl  •t-lw2<E/I  ■VAX.lJX.CX,UX.AY.UY.CY.DY.AZ.b2«C2.D2/2/Kl 

4  /6/NXl  •NX.'  •NX3«NY  1  .NY2  4NY3«N2l4N22»N23 
real  LL  •'_•<! 

real  NXl •NX2.NXj4NYi •NYa^NYO^NZl IN224NZ3 
READ  994N 
00  9  I  I»1 .N 

50  read  1 104AnaME4UL4LL4L4DL4DU4DD4AX4DX4CX4DX4AY40Y4CY4DY«AZ4DZ»CZ. 
1  02 

OELK3&.2e3ia/L 

E-DELK/50. 

OE'  <K»2.  jE-5/L 
I  2  I  »  1  OU-OL  1  /OO*  I  • 

DO  15  I.?2sl4t21 
DOO=OL+i:U*(  lZZ-1  1 
PRINT  1C04ANAME4UL4LL4L 
PRINT  1024000 
PRINT  101  •E^Ob.LK^OELKK 
00  1  1=1-2 
W I >6943. 

IFtl.EQ.2>  Wl=10600« 

NXl »FIN(W) 4 11 
NY1»F|N(w1 .21 
N2i=FlN(iKl4  3l 
W2«LL 
GO  TO  11 

10  w2*w2-w2/8» 

SF(w2.1.T.ULI  go  TO  1 

11  W3«  I  •/(  1  4/11(1  + 1  ./WZI 
IF< wl 4GT4LL.OR.Wl .LT.UL.OR.w2.GT. 


LL.OR.W2.LT .UL . OR . W3 t GT . LL . OR . W3 


00000800 
OOOOObOO 
00000700 
00000800 
00000900 
00001000 
00001 100 
00001200 
00001300 
0U001400 
00001500 
00001600 
00001700 
0000 1800 
00001900 
00002000 
00002100 
00002200 
00002300 
00002400 
00002500 
00002600 
00002700 
00002800 
00002900 
00003000 
00003100 
00003200 
00003300 
00003400 


12 


.lt.uli  go  to  10 

00003500 

NX2«FIN(w24 1 1 

00003600 

NY2*FIN(WE421 

00003700 

NZ2«F INIW2431 

■-10003800 

NX3=F IN( W3. 1 1 

.10003900 

NY3«FIN(W3421 

00004000 

NZ3=FINIW3.3I 

00004100 

0IV»D0D/I  (NXl+NZI  1.'2. 1 

00004200 

0«ASIN( (W2/w3l*5IN(0IV  ♦ASIN( (W3/W2l*SINI-OIV  )1I1*C1 

00004300 

DO  4  K> 1 .6 

00004400 

1F(K.GT.3|  GO  TO  6 

00004500 

In3s2 

00004600 

INI>1 

00004700 

INZal 

00004800 

1F(K.E0.2)  IN1«2 

00004900 

1F(K.E0.3I  lri2»2 

00005000 

GO  TO  12 

OOOCSIOO 

1N3.  1 

00005200 

INI=2 

00005JOO 

1N2  =  2 

00005400 

IFik.EO.SI  IN1>I 

00005500 

IF(K.E0.6>  IN2>1 

00005600 

Jr-l 

00005700 

TCL«FNX<0..0. .0. 1 

00005800 

IF( IX.GT.l i  GO  TO  5 

00005900 
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I 

IS 
9 
99 
1  00 


CALL  HwOK(U*TCL«0.«4) 
ir( IXtST. 1 )  GO  TO  S 
HWCl-HWlACl 

Hwca«Mw2»ci 

FN»FN2<  TCL» 

FXHW1«AS1N(SIN(HW1 )*FN)*CJ 
EXHwSa  AS  I N( S  t  N I HW2  >  AFN ) »C 1 
PClaPVAR(TCU<0«>0tV/2«  <0* ) 

PC2«PVAR(TCI..0.  .-OIV/2*  lO.  > 

CALL  HWDX»LtTCL+DIV/2..PCI «6t 
hkC3«Mw2 

CALL  HwDK<L»TCL-0IV/'2*  .PC2«5) 

HWCA»MW I 

A 1  -0  I V/2  •  ♦A6S  ( PC  I  -f  MWC3 ) 

A2«0  I  V/2  •  «AeS  (  PC2'fHWC4  ) 

FN«Fn2(TCL) 

FXWl ■ABSt ASIN( SItJ(  Al  )*FN)+ASIN«SIN<  A2  )*FN))*C1 

MWC3«HWC3*C1 

MWC<»»HWC4*C1 

FNH«<AJ*A2)*CI/'0  - 

OQ-O* 

TNCL»  TNCPM ( TCl  »  00 ) 

CALL  HW0K(L« TNCL.PNC<3) 

HWMC1«HX1*C1 
HWNC2«HW2«C1 
FN«FN2 ( TNCL+PNC 1 
eXMWS«ASlN<SIN(HM(l  )*FN»*Cl 
exMwe* AS I N( s I N<  >  *fn  » *c i 

)F(j.ea.2i  GO  TO  e 
CALL  OUT 
J>2 

00»DEL< 

TCL«FNX<0«<0<  <0Q) 

GO  TO  7 

PNC  1  ap VAR '!  TNCL  «PNC  <  0  I  V  tOELK  ) 
call  HVOKCL* TNCL+OI ViPNCJ .3) 

HWNC3«HWI 
hwnca«hik2 
FN»FN2 ( TNCl+PNC ) 

EXW2»ABS<*SIN(SIN»MWNC3l*FN»-ASlN(SINIMWNCA )*FN) >*C1 

liWNC3aMWNC3*Ct 

MWNCA«MWNCA*Cl 

FNL2»AeS(HWNC3-HWNCA  > /D 

CALL  Out 

continue 

GO  TO  10 

CONTINUE 

continue 

CONTINUE 
FORMAT! 12) 

FORMAT!  20X. 2 IM***  UP-CONVERSION  IN  tAStAH  *»*^/39M  THE 


00006000 

00006100 

000061200 

00006300 

00006400 

00006500 

00006600 

00006700 

00006600 

00006900 

00007000 

00007100 

00007200 

00007300 

00007400 

00007500 

00007600 

00007700 

. .  00007600 

00007900 
00008000 
00008100 
00008200 
00008300 
00008400 
00008500 
00008600 
00008700 
00008800 
00008900 
00009000 
00009100 
00009200 
00009300 
00009400 
00009500 
00009600 
00009700 
00009800 
00009900 
oooioooo 
00010100 
00010200 
00010300 
000 10400 
0001 0500 
00010600 
00010700 
00010600 
00010900 
SHORTEST  WAOOOHOOO 


IVELENOTH  transmitted  IS  «F10*223eM  THE  LONGEST  WAVELENGTH  TRANSM I TOOOl 1 1 00 
2TED  IS  «F10.2/2SM  THE  LENGTH  OF  CRYSTAL  IS  .F 1 0 *2* 2X t 3HCM./  90H  AL00011200 
3L  PHASE  MATCH  ANGLES  ARE  IN  OEGREES«  P»M.  ANGLE.S  ARE  MEASURED  8ETW00011300 
4EEN  THE  Z-AXIS  ANO  Kl  /’33H  TCL  ■  COLINEAR  PHASE  MATCH  ANGLE  /'37M  T00011400 
SNCL  ■  NONCRITICAL  PHASE  MATCH  ANCLE  /-AOH  PNC  ■  ANGLE  BETWEEN  i  Kl  tKOOOl  1500 
621  FOR  NONCRITICAL  P*M.  /36H  ALL  HALFWIDTHS  ARE  IN  MILLIRADIANS  /  00011600 
777M  HwC  «  MALFWIDTH  OF  DELK  VS.  THETA  CuRVE  MEASURED  BETWEEN  MAX.  00011700 
band  first  2ER0/39M  HWNC  »  HALFwIOTH  FOR  NONCRITICAL  P.M.  )  00011800 
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102  FORMAT*  68H  PCI.PC2  *  PCI  FOP  TCL  »  TCL+OR-O I V/2 .  TMC  INTERNAL  UtVOOOMVOO 
lERGF.NCE  OF  W1  /  31H  PNC  J  •  PNC  FOP  TNCL  ■  TNCL+DIV  /■  32H  Ew  •  EXTE00012000 
JRNAL  ANGULAR  APERATURE  /  32H  R  >  NUMBER  OF  RESOLVEABLE  LINES  ✓  00012100 

349H  EXTEHNAL  divergence  (FULL  ANGLE!  OF  PUMP  BEAM  »  .F6.A)  00012200 

101  FORMAt(46H  MWC  ANO  HWNC  are  measured  with  TCL/TNCL  fixed  /l^H  UELT00012300 
lA  K  IS  WITHIN  .F7.3,13H  IXCM.  OF  2PI/'L/'<SH  2PI/L  *  .F7.3.SX.2AH  <2P000I2400 
2I/’L)/K3  IS  APPROX  K  •E10*3/46h  K1iK2«  ANO  K3  ARE  COPLANEP  IN  THE  (00012500 
3T«Z)  PLANE  /37H  E  POLARIZATION  IS  IN  THE  (V<Z>  PLANE  <'41H  O  POLAR 1000 12600 
4ZA7I0N  IS  OUT  OF  T(4E  (Y«Z)  PLANE  /69H  FOR  UNIAXIAL  CRYSTALS  THE  Z-QOO 12700 
SAX  IS  ■  C-AXIS<  NZ  >  NE>  ANO  NX  •  NY  «  NO  xESH  •  INDICATES  DELTA  K  00012000 
6>i  ZPI/L  /-SSH  wl<W2<W3  correspond  TO  PUMP,  IR.  AND  SIGNAL  WAVELENGTOOO 1 2900 
7HS  /SSH  NX,NY,N2  ARE  REFRACTIVE  INDICES  FOR  X,Y.Z  OPTICAL  AXIS  //500013000 

fl9H  ALL  POSSIBLE  uP-CONvERS  I  ON  PROCESSES  ARE  LISTED  BELOW - ✓/’/’/)  000 1 3 1 00 

110  FORMAT(A5,5X,6F10.1/4F10.4/'4F10.4/4F10.41  00013200 

END  00013300 

function  TNCPM(TCL,0)  00013400 


COMCdON  ✓lO/'P 

T»TCL 

OP=O.05 

5  T1 -FNX(DP. 3. 14159.0) 

IFITI.LT.T)  go  to  10 
OP«-OP 

T1 ■rNX(DP»3. 14159,0) 

IF(Tl.LT.T)  GO  TO  10 
0P«DP/10. 

|F(ABS(0P>,LT. 0,00001 >  CO  TO  21 
GO  TO  5 
10  ’■•Tl 
PbOP 

14  Tl  ■FNX(P-*-0P,3»  14159.0) 
IFITI.LT.T)  GO  TO  12 
p.pWOP 
OPa-OP/4, 

1F(ABS(OP>,LT.O.OOOI )  CO  TO  20 
T«T  I 

GO  TO  14 
12  P«P+DP 
T«T1 

GO  TO  14 

20  TNCPM.TI 
RETURN 

21  TNCPM.TI 
P»OP 
RETURN 
END 

FUNCTION  PVAR( T.P.OD.E ) 

COMMON/ 1 /w 1 , W2 , W3/2/X I 

real  KI 

D*^  .001 

1F<OO.GT.O.O)  0—0 


00013500 
000 1 3600 
00013700 
00013600 
00013900 
00014000 
00014100 
00014200 
00014300 
00014400 
00014500 
00014600 
00014700 
00014800 
0001 4900 
00015000 
00015100 
00015200 
)00 15300 
00015*00 
00016500 
00015600 
00015700 
00015800 
00015900 
00016000 
00016100 
00016200 
00016300 
00016400 
00016500 
00016600 
00016700 
00016800 


PP*P 
Tl  bT+OO 

<1 ■6,2831 eE8*FNI ( Tl )/Wl 

1  PP  =  PP+D 

2  X1=ABS(FDK(TI .PPITE) 
x2-ABS(FDK( T 1 ,PP+D1+E1 
if(x2,lt.xi)  go  to  1 
ppepp.D 

D«-D/5.0 


00016900 

00017000 

00017100 

00017200 

00017300 

00017400 

00017500 

00017600 

00017700 
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IF( ABS(O) .LT. lE-5)  GO  TO  10 
GO  TO  2 
10  PVAR=PP 

return 

END 

FUNCTION  FNX«XX»PP.0D) 

data  (E»l4E-6) 

COMMON/ 1 /■W]  .W2,W3  /4/JJX 

Ft<(0,A»=SO»T( <FNt<A)/wl )**2+FN2C«Q)/(W2*W2)-2.*FN1 (A)* 
1  SORT (FN2C I  0) > *CXX/ ( Wl •W2) ) 

0=OD/6.2a3iaEB 

P=PP 

x=xx 

IF(XX<GT.O)  GO  TO  20 
X»-X 

P*  "J  .  I  A  I  !S927«  +  PP 
20  IKX=|Jx=0 
01 »AeS( XP-XX ) 

XP  =  XX 
CX=COSlX> 

SXeSlNlX) 

CXX=C0S(3. 1A1S927-X) 

CP=COS<P> 

^ ^ * R* E T »0 407# OR« R4GE *  1 «S7 )  R*0407 
IF  I O I 4LT4 1 00 4 ♦£ •  01>1004*E 
IF (01 4GT4O4 I  I  01«04l 
IF <P4EQ40 407 )  OICO4I 
1  IXX=|XX+1 

1F( IKX4GT4 100 )  GO  TO  5 

R0*R*Ol 

CR*COS(R) 

SR«SIN(P) 

CR0»C0S(RO) 

SR0=SlN(R0) 

C2i!CR*CX+SR*SX*CP 

C5«CRO*CX*SRD*SX*CP 

YI=FX(C24R) 

V2=FK<C54R0I 

C3«IFN1 (Rl/Wl-CxX*SORT (Fn£C(C2> )/w2)/YI 
CA»(FN1 (RO)/Wl-CXX«SQRT(FN2CtC5> )/W2|/V2 
C6=C3*CR+SR* sort I ABS  <  1 4 -C3*C3 »  »  *CP 
C7*CA*CRD+SRO*SQRT ( A0S  < I 4-CA*CA  >  »  *CP 
Y» ABS  <  SORT  <  FN JC I C6  > I /w3T0”Y 1 ) 

2=ABS< SORT ( Fn3C ( C7) I/W3+0-Y2) 

IF  (A0S(DI)4LT,  e  1  GO  TO  2 
iF( ixx4Ea4 1 )  GO  TO  10 
IF  <Y4GT.2)  GO  TO  3 
01 =-Dl/l0 
R*R-10*0I 
GO  TO  1 
3  R>R+DI 

■  F<R4l.T.04040R4R.Ge4  1  >6  )  Go  Tu  12 
GO  TO  I 
2  FnxcR 

IF(C64LT.04)  IJX«2 
RETURN 

10  IF(Z4GT4Y)  0I=-D1 
IF<Z.GT4YI  GO  TO  1 
GO  TO  3 


00017600 
00017900 
00018000 
oooiaioo 
00018200 
00018300 
00018AOO 
000  1 8500 
00018600 
00018700 
00018800 
000 18900 
00019000 
00019100 

00019200 

00019300 
00019400 
00019500 
00019600 
000 19700 
000198C0 
00019900 
00020000 
00020100 
70020200 
00020300 
OOOaqAoO 
00020800 
00020600 
00020700 
00020800 
0002o9o0 
OOO2100O 
00021100 
00021200 
00021^00 
00021400 
00021500 
00021600 
00021700 
00021800 
00021900 
00028000 
00022100 
OOO222OO 
00022300 
00022400 
00022500 
00022600 
00022700 
00042800 
00028 vOO 
00023000 
00023100 
00023200 
00023300 
00023400 
00023500 
00023600 


I 


» 
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12  tF(iJX«aT«0)  00  TO  6 
14  IJX>1 
RaO.07 
OlaO* 1 
IXXal 
GO  TO  I 
6  IJXa2 
RaO.O 
return 

9  IF< I JX.GT.O)  00  TO  e 
GO  TO  14 

S  PRINT  1444X<P«R 
I  JXa2 
PNXaR 

144  FORMAT* lOH  FNX  ERROR  .5X«3HX  a  4F9.a«bX«3HP  *  .F944.iX.3HR 

RETURN 
END 

FUNCTION  FOX*  T.P) 

COMMON  /1/Wl  .W2<V*3/2/Kl 
REAL  K2»K2S«X2C«X3.K| 
data  (Cl>6.28318Ee) 

K2aCl*FN2<T+P )/W2 
K2SaK2»SIN(P) 

K2CaK24C0S<P) 

P3a ATAN ( K2S/ <  K 1 +K2C >  > 

K3aC  1  apNO  (  T+P3  )  /IK3 
FDXaK3-saRT<K2S**2+*X2C+Kl ) **2 ) 

RETURN 

END 

SUeROUTINE  HV*OK<L.T«P.  I  ) 

real  L«L2«Kt 

DATA  <CCa«>2e3l6Ee) 

COMMON  /1/XI.W2.W3  /2/'Xl/5/Vl  .V2»E/«/lX 

KlaCC*FNl «T»/Wl 
L2a  64283ie/L 

IF(ABS(F0K(T«Pn.GT4L2*l.S  )  GO  TO  30 
OVaOaOl 

IF* I .EQ.E.OR. I .EQ.6I  OVa-.Ot 

1  VaOV 

2  OKaFOX* T.P+V) 

IF(I.GT.3)  OKaABS(OX) 

IF<At>S<  OX  -L2)4LT4EI  go  TO  20 
IF*  OX  .LT.L2>  GO  TO  10 
DVaOV^lOi 

IF(ABS*0V1.LT.IE-6)  go  TO  31 
VaV-0V*9. 

GO  TO  2 

10  VaV+OV 
GO  TO  2 

20  IFIDV.LTtO. >  GO  TO  29 
OVa-O.Ol 
VI  av 

IF ( 1 .EOo 1 .OR. I .EQ.S)  GO  TO  25 
GO  TO  1 

29  V2av 
RETURN 

30  IXa2 

31  PRINT  100.T.P.L2 
100  F0RMAT(3F1<!.SJ 


00043700 
00043BOO 
0002.  >00 
0002-.000 
00024100 
00044200 
00024300 
U0044400 
00024500 
00044O00 
00024700 
00024800 
00024900 
00025000 
= .F9.4 ) 00045100 
00045400 
00025300 
00045400 
00023500 
00025600 
00025700 
00025800 
00045900 
00026000 
00026100 
00026200 
00026300 
00026400 
00026500 
00046600 
00026700 
00026800 
00026900 
00027000 
00027100 
00047200 
00027300 
00027400 
00047500 
00047600 
00047700 
00047600 
00027900 
00028000 
00028100 
00028200 
00028300 
00028400 
00028500 
00026800 
00028700 
00028800 
OO0269O0 
U0029000 
00029100 
UU029400 
00029300 
00029400 
00029500 
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RETURN 

END 

subroutine  out 

common/ 1 /"wi  .WE tw3  /4/1x/:o/pn 

1/1  I /K . TC. TNC .HW 1  .HW2. MWNl  ,MWNE. hw3 . HW4  »MWN3 . HWN* .UU . J  «UW 
2/lE/EMWl  .EHWE.EXWl  .EXW2.EHWb.EHW6.KNLl  .ENL2.HLL1  .HNLLl  .t->LL2 
3  /6/Nxl .NX2.NX3.NVI ,Mv2.NY3.NZI .NZ2.NZ3 
real  NXI .NX2.NX3.NV1 .NV2.NY3.NZ1 .NZ2.N23 
DATA  <uP«57,29B77951 . (ON’. 01 7453293) 

IF ( K.GT. 1 .OR. J.GT. 1 >  GO  TO  5 
01 V’OV*1000. 

I  PRINT  lOO.Wl .W2.W3.NX1 .NX2.NX3.NY1 .NY2.NY3.NZ1 .NZc.NZJ.UU.uI V 
5  IF ( IX. GT. 1 )  GO  TO  10 
TCL»TC»UP 
PC1=PCLI*UP 
PC2=PCL2*UP 
TNCL»TNC*UP 
PNC»PN*UP 
PNC  1 «PNCL 1 *UP 
IF(J.E0.2)  GO  TO  3 
( 


. 12. 13.14.15. 16) .K 
110. TCL .MW  1  .MW2 . EMW 1  .EHW2 . TNC L . PNC . MWN 1 , HWN2 . EMW5 .EMW6 

120.TCL.HW1 .HW2.EHW1 .EMW2.TNCL.PNC.HWN1 ,hwN2,EHW5.EHw6 

130.TCL.HW1 .HW2.EHW1 .EKW2.TNCL.PNC.HWNl . MwN2 , EHW5 .EHWb 
4 

140.TCL.HW1  .MWZ.EHWl  .EHW2.TNCL.PNC.HWN1  . HWN2 . EHWt>  . LHW6 
t 

150.TCL.MWI .MW2.EMW1 .EHW2.TNCL.PNC.HWN1 .HWNZ.EMWb.tMWO 
I 

150.TCL.HW1 .MW2.EHW1 . Ehw2 . TNCL .PNC . MWN 1 , MWN2 .EMW5 .EHW6 


GO  TO 

11  print 

RETURN 

12  PRINT 
RETURN 

13  PRINT 
RETURN 

14  PRINT 
RETURN 

15  PRINT 
RETURN 

16  PRINT 
RETURN 

10  GO  TO  (21 .22.23.24.25.26) .K 

21  PRINT  210 
RE  TURN 

22  PRINT  220 
RETURN 

23  PRINT  230 
RETURN 

24  PRINT  240 
RETURN 

25  PRINT  250 
return 

26  PRINT  260 
RETURN 

3  PRINT  170.  TNCL.PNC.HWN1 .HWN2,tHW5.EMw6 

PRINT  200. PCI  .PC2.HW3.HW4.EXWl  .  FNL  I  .PNLl  .  HVt'N3  .  HWN4  .  EXW2  .FNL.2 

if(k.eo.5)  go  to  4 

RETURN 

4  PRINT  300 
RETURN 

100  F0RMAT(6H  Wl  =  . F 1 O . 2 .  1  OX . 5HW2  »  . F 1 O. 2 .  1 CX . 5Hw3  ’  .F10.2.I0X/3 
I  NX  =  .F10.4.9X)/J(6H  NY  «  .FI0.4.9X)/3(6H  NZ  «  .F10.4.9X)// 

25X.  35HINTeRNAL  RESOLUTION  LIMIT  (MRAO)  =  .F7.4. 


3  5X.  16H1NTERNAL  OlV. 

110  FORMAT (9M  0+0=E - .7M 

1  F7.4.2X.6HPNC  = 

120  format (9M  E+O’E - .7h 


’  .F6.4//) 

TCL  «  .  F7.4 .2X.6HHWC  « 

.  F7,4 .2X.7HHWNC  ’  .4F6. 
TCL  «  •  F7.4.2X.6HHWC  = 


.  4F6,  I .2X.7HTNCL  » 
I  I 

i4F6. 1 .2X.7HTNCL  ’ 


00049600 
00049700 
00049600 
00049900 
00030000 
00030100 
00030200 
00030300 
00030AOO 
00030POO 
UOOJ0600 
00030700 
00030B00 
00030900 
00031000 
00031 100 
00031200 
00031300 
00031400 
00031500 
00031600 
00031700 
00031800 
00031900 
00032000 
00032100 
00032200 
00032300 
00032400 
00032500 
00032600 
00032700 
00032800 
00032900 
00033000 
0U033100 
00033200 
00033300 
00033400 
00033500 
00033600 
00033700 
00033800 
00033900 
U0034000 
00034100 
OOOJ4200 
00034JOO 
00034400 
00034500 
00034600 
00034700 
(  6HOOO->4eOO 
00034900 
00035000 
00035100 
»  .00035200 
00035300 
.00035400 
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I  F7.4i2X.6MPNC  •  F  4.CX.7HHWNC  ■  .<il-6.l) 

130  rORMAKVH  07E«t.~— «7m  TCL  •  f  7  ,  A  1 2X*6MHwC  ■  «  AF6.  I  «  ix  .  7(iTNtL  • 
\  F7,*,2X.6MPNC  •  F7. 4  tix  .  /i-.«WNC  ■  .AF5.  I  ) 

140  FORMATION  E+E^O - -  7H  Tex  s  <  F7  A  2X  '  tiHHMC  ■  •  AFfc ,  1  .  .  7HTNCL  ' 

I  F7.4*2X.6MPNC  ■  «  F 7.rt ,2x , 7MMWMC  •  .AF6.1) 

190  formation  0-fE«0- — «7H  TCI-  «  .  F7.  •.  ,2x,6H  IwC  »  .  AF6.  1  .  2X  .  7MTNCI-  = 
1  F7.4,2X*6MPNC  ■  I  F7.A.2x.7MNyiNC  «  .aF6>I) 

160  formation  EA0>0 — -t7M  TCL  ■  «  F7«A,2X  ©nhwC  •  .  AF&.  1  •  <dX .  TriTNCi.  * 
I  F7«4«2X<6HPNC  >  •  F7.4.3Xt7F-lwNC  ‘  «AF6«n 

170  FORMATION  •  <  7HTNCL  = 

I  FT.*»2^{.  6riKNC  •  F7,  T»2X«7NNWNC  «  .AF6.  I  ) 

200  FORKA' < 10X«SNPC1  ■•F7.4,6H  PC2  •«F  ,4«bH  NW  «.2F6.\.AH  tW*.F6.l. 
1  3N  Hr,F5.'. -7H  PNCI  •«F7,A(ftH  Nw  k«2K«,1.AN  EW«.F6.1f3N  P:.F6. 


Phase.  mATChAole  **« 
P'  ASfc  MATCNAcslE  **« 
P  .ASfc  MATCHAbLt  *F* 


210  /•■0RM«r(9H  0+0«E-'-*  not  PHASt  MATCNAciLt  *** 

220  FORMAT(ON  £>0*8 - •10Xi2eH»*»  NOT  PHASE.  MATChAblL  *** 

230  FORMATION  O+E-E — -«10X*2eH*«*  NOT  PHhSE-  MATChablc  *** 

240  formation  E+E«0 — -•10X«2eH**F  NOT 

250  FORMATION  0+E-O— .  1  OX  •20N*»«  NOT 

260  formation  e*0«c - not 

1  //) 

30'!  FORMAT! //aoi  3M  *!✓/» 

END 

FUF.CTlON  F1NHK»  1  » 

COMMON/ 1  3/AX  (BXiCXtOXt  AT  «ey<CY«Oi'<AZ,B2«CZ.OZ 
X»W*W*l.OE-e 
GO  TO  ( 1 <243) <  t 
J  FIN»SQRT( AX*ex/(X-CX)-OX*X 
return 

2  F1N«S0RT( AY+BY/(X-CY)-DV*X) 

return 

3  FlN«SQP<  <A2*t#2/(X-CZ> -rJIAXl 
return 
END 

FUNCTION  FNl«A) 

common  /1/W14W24W3  /O/lNl 

4  /6/NX  I 4NX24NX34NY I 4NY2  4NY34NZ 1 4NZ24NZ3 

real  NXl tNX24NX34NYl 4NY24NY34NZ1 4NZ24NZ3 
GO  TO  ( 1 .2) 4  INI 

1  FNl-NXl 
RETURN 

2  FNI«1»/SQRT(  ICOSI AI/NYI )**2+(S1N«a>/nZ1 >**2) 

RETURN 
END 

FUNCTION  FN2IA) 

COMMON  /1/W14W24W3/8/IN2 
4  /6/NXl 4NX24NX34NYI 4NY2 4NY34NZ 1 4NZ2 ,n23 
real  NXl 4NX24NX34NY1 4NY24NY34NZ1 4NZ24NZ3 

GO  TO  (1 42i 4 1N2 

1  FN2aNX2 

return 

2  FN2«1 ./SORT! !C0S( A)/NY2)**2F(SIN|A)/NZ2»**2> 

RETURN 
ENO 

FUNCTION  FN3(A) 

COMMON  /1/Wl 4W24W3/7/IN3 
4  /6/NXl 4NX24NX34NYI 4NYi:  4NYJ4NZ 1 4NZ2,nZ3 
real  NXl 4NX24NX34NYI 4NY24Ny3,NZl 4NZ24NZ3 
GO  TO  ( 1 42) 4 IN3 
1  FN3«NX3 


) 

) 

) 

) 

) 

//AVON 


OOOASbOC 
4CC0-'bPOC 
0003S7CC 
4  OOCAbbOO 
OOOA5900 

4  D00-J6000 

000-56100 
4C'OC-56cOO 
C00-56i00 
4  000-56400 
000A6b0C 
000A6600 
1 )OOOA6700 
000->6b00 
UOO-56VOO 
000-57000 
OOOA7100 
0O0-57C0O 
♦ )OOOA7AOO 
C00A7A00 
00037500 
000-57600 
000.57700 
OOOA7aOO 
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RETV^N 

2  FN3*1  ./.5QRT(  (COS(  A)/'Ny3)**2*(SIN(A)/NZ3l**2l 
RETURN 
END 

FUNCTION  Fn2C(Q) 

COMMON  /l/Wl «W2»W3/B/|N2 
A  /6/NX\ .NX2»NX3«NY1 .NV2«NY3«NZ1 .NZ2.NZ3 
REAL  NXl «NX2«NX3*NY1 .NY2.NY3.NZ1 .NZ2,NZ3 
GO  TO  <1 <2) « IN2 

1  FN2C»NX2*NX2 
RETURN 

2  FN2C* I ./( (Q/NY2 1**2+ ( 1 .-Q*« I /NZ2**2 ) 

RETURN 

ENO 

FUNCTION  FN3C(0) 

COMMON  /l/wl  ,w2«\M3/-7/IN3 
A  yG/Nxl <NX2«NX3«NY1 «Ny2>NY3*NZ1 «NZ2»NZ3 
REAL  NXl <NX2»NX3«NY1 «Ny2»NY3«N21 iN22,N23 
GO  TO  ( I iZ) < IN3 

1  FN3C*Nx3*Nx3 
RETURN 

2  FN3C«1 .Xl (0/Ny3)**2+< l.-a*0*/N23**2) 

RE-URN 

ENO 
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